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We study the rho meson mass in a uniform background magnetic field eB at zero temperature, in 
search of indications for the magnetically induced rho meson condensation, as predicted recently 
by Chernodub. The holographic model used is the Sakai-Sugimoto model with two flavours 
and a non-zero constituent quark mass. We fix the free holographic parameters by matching 
them to the phenomenological value for the constituent quark mass and the experimental values 
for the pion decay constant and the rho meson mass, this in absence of a magnetic field. In a 
first approximation, the Landau levels are recovered, indeed indicating an instability of the QCD 
vacuum at a critical magnetic field, eBc ^ m^p,Xo a. phase where rho mesons are condensed. We 
improve on this result by also taking into account the holographic analogue of chiral magnetic 
catalysis, numerically solving the mass eigenvalue equation for the rho meson, which depends 
on eB both explicitly and implicitly through the changed embedding of the flavour probe branes. 
This turns out to raise eBc with a few percents. As a byproduct of our analysis we find that the 
separation between the chiral symmetry restoration temperature T^{eB) and the deconfinement 
temperature is 3.2 percent at eB = 30m^ w 0.57 GeV^. 
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1. Introduction 

The interest in QCD in the presence of a strong magnetic background has arisen recently since 
strong magnetic fields are expected to appear in heavy ion collisions [1]. In [2, 3] it was suggested 
that the cold {T = 0) QCD vacuum is unstable towards formation of a condensate of charged rho 
mesons when a sufficiently strong external magnetic field' (S ~ wip) is applied. The argumentation 
in [2] is based on the DSGS-Lagrangian [4] describing self-consistent quantum electrodynamics 
for the rho mesons, whereas in [3] a similar- reasoning was built up in a NJL-model. 

Inspired by these papers and the already known existence of rho meson condensation at finite 
isospin chemical potential in holographic QCD [5], we have studied the dependence of the rho 
meson mass on a background magnetic field in a holographic approach. Holographic QCD-models 
give a description of hadronic physics through a dual supergravity theory in a higher-dimensional 
world. The duality is vaUd at large 't Hooft coupling where QCD itself is unmanageable, thus 
providing a setting for studying non-perturbative QCD effects. The holographic model which we 
have used, the Sakai-Sugimoto model [6, 7], manages to reproduce much of the low-energy physics 
of QCD, in particular confinement and dynamical chiral symmetry breaking at low temperature. 

In section 2 we give a short review of this model, including our numerical values for the free 
parameters obtained by matching to phenomenological and experimental results, and we describe 
the method to introduce the magnetic field. The effect of the magnetic field on the geometry of 
the model in the current approximation is summarized in section 3, and then used in section 4 to 
determine the B-dependence of the rho meson mass and the critical value of the magnetic field at 
which the charged rho mesons condense. In the last section we calculate the magnetically induced 
split between deconfinement and chiral symmetry restoration temperature, finding a reasonable 
agreement with phenomenological and lattice output. In this proceeding, we only summarize our 
findings, details will appear shortly in [8]. 

2. Holographic setup 

The Sakai-Sugimoto model [6, 7] involves a system of Nf pairs of D8-D8 flavour probe branes 
in the D4-brane background 

ds^ = ['^f\il,ydxV+f{u)dr')+(^^y' + , 

• ^'-t'- ^w = '-i' 

where dQ.^, £4 and V4 = 871^/3 ai^e, respectively, the line element, the volume form and the volume 
of a unit four-sphere, while /? is a constant parameter related to the string coupling constant g^, the 
number of colours Nc and the string length 4 through = ngsNci^^. This background has a natural 
cut-off a.tu = UK and the QCD-like theory is said to "live" at m — 00. We set Nc = 3, and the number 
of flavours Nf = 2. We stress here that we ignore the back reaction of the flavour branes on the 
background, meaning we are working in a holographic analogue of the quenched approximation. 

'For notational convenience, we shall usually write just B, rather than eB. 
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Figure 1: The Sakai-Sugimoto model. 



The Sakai-Sugimoto model incorporates dynamical chiral symmetry breaking as a conse- 
quence of the ?7-shaped embedding of the flavour branes: the D8-branes and D8-branes merge 
at a certain value of the extra holographic dimension u = uq. In contradistinction to the original 
setup of [6, 7], we consider the more general setting with uq > uk, this in order to have non-zero 
constituent quark mass [9] 

1 /•"o du 



Ina' 



(2.1) 



with a' = the string tension. The bare quark masses are always zero, so we are working in the 
chiral limit. 

The gauge field A,„{x^,u) (m = 0, 1,2,3,m) living on the D8-branes describes mesons in the 
boundary field theory. The action for this gauge field is given by the non-Abelian DBI-action^ 

Sdbi = -TsJ d^^dx £4e-^STr^ - det [gZ + (27ra')^m„] (2.2) 

with T% = l/{{2n)^£'^) the D8-brane tension, STr the symmetrized trace, g^^ the induced metric 
on the D8-branes, and F,„„ = d,„An — dnA,„ + [A,„,A„] = F^f^j" the field strength with anti-hermitian 
generators 

f^' = ^(1,01,02,03), Tr(?V) = -^. 
The parameters R,gs,is,MK, uk and A = gyM^c are related through the following equations: 

^' = 0^^ S^ = ^4f^ ^K = IWI (2.3) 
2Mk 2kMkIs 9 

Without loss of generality one can put uk= ^ [7]. To fix the seven unknown parameters (/?, A, 
Is, Mk, Uk, gs and uq) we thus need three additional conditions. We intend to choose these in such 
a way that the major features of Nf = 2 QCD are reproduced as well as possible. In particular, we 
use the phenomenological value for the constituent quark mass {niq = 0.310 GeV, as in [10]), and 
the experimental values for the pion decay constant (/^ = 0.093 GeV) and for the rho meson mass 



^Because the 't Hooft coupling X is assumedly large, we ignore in a first approximation the Chem-Simons part of 
the action in the analysis, being a factor 1/A smaller than the DBI-action. 
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in absence of magnetic field {nip = GeV). The results of our numerical analysis are [8] 

Mk ~ 0.721 GeV, mq ~ 1-92 GeV^^ and K = ^ 0.006778. (2.4) 

From this, we do extract a relatively large 't Hooft coupling, A ^ 15. Similar values can be found 
in [7], however these were obtained in the limiting case uq = uk, corresponding to zero constituent 
quark mass. In this setting it would be, as we shall see shortly, impossible to model a magnetic 
catalysis effect, which is by now generally accepted to occur in QCD. 

To turn on an electromagnetic background field A^'" in the boundary field theory we use the 
standard technique [7] of putting 

A, (« ^ oo) = = el^'^ ^ A- = e + las) (2.5) 

which, for the case of an external magnetic field in the ;C3 -direction in the boundary field theory 
{F^f = 5jA|'" = B), amounts to setting 

A^{u oo) = e2e„,xiB5^2 = ^^^^^^^ ("^) +-^'^^^^^2 (^"y") ' 
or, adapting the notation A^ = A^ (m — )• oo) : 

A\=xieB and A^=A2/3. (2.7) 
Because we want to look at the effect of the electromagnetic field on the rho mesons, which trans- 

3 

form in the adjoint of U {Nf)i x U {Nf)R and hence only couple to the isospin component A^ of the 
electromagnetic field, we may ignore A^ . 

We work in the A„ = gauge and fix the residual gauge symmetry by working in a particular 
gauge where the total gauge field is given by [7] 

A^ (x^ , „) = A^ + £ {x^' )Yn{u), (2. 8) 



with V^,,i{x^) a tower of vector mesons with masses and {Yn{u)}n>i a complete set of functions 
of u, satisfying the eigenvalue equation 

«'^'7e '^'(")^« [u'^^rB'^\u)duYn{u)] = -R^mlYniu), (2.9) 
and the normalization condition 

1 / 1 

— {2Ka')2R^^^TsV4 / du u-^/^Yl^^^R\,„{u)Yn{u) = -d,nn, (2.10) 

Jun 4 



2g. 
with 

rB{u)= , \, r ^ A m =1+b' - , Ao = l+52 - . (2.11) 
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3. Influence of the magnetic field on the embedding of the flavour probe branes: 
magnetic catalysis of chiral symmetry breaking 

As already discussed in [11] for the Sakai-Sugimoto model, an external magnetic field pro- 
motes chiral symmetry breaking: when holding the asymptotic sepai^ation L between D8- and D8- 
branes fixed, uq, which is in a one-to-one correspondence with L [6, 7, 9], grows with B. This 
means that the probe branes get more and more bent towards each other, driving them further and 
further away from the chirally invariant situation (straight branes). This feature corresponds to a 
holographic modeling of the magnetic catalysis of chiral symmetry breaking [12]. For our value of 
L = 1.57 GeV^^ (coiTcsponding to uq{B = 0) = 1.92 GeV^^) kept fixed, the numerically obtained 
B-dependence of uq and consequently of the constituent quark mass m^, via the relation (2.1), is 
plotted in Figure 3 [8]. Since niq rises with B, we expect that taking this magnetic chiral catalysis 
into account will translate into the rho meson mass mp also rising with B, at least when ignoring 
the lowest Landau level shift (see next section). 



Mo (GeV"') (GeV) 




Figure 2: mq as a function of the magnetic field. Figure 3: The constituent quark mass as a func- 
tion of the magnetic field. 



4. Equations of motion for the mesons in the magnetic field 

Plugging the gauge field ansatz (2.8) into the action (2.2) and expanding it to order V^,,, we 
derive the equations of motion for the vector mesons „, and thus in particular the rho mesons 
(we concentrate on the lightest vector mesons, because they will probably be the first to condense 
at strong B). We find [8] 

mlp^-Dlp^+HG^vP^- =0, withD^=5^-/M;;'", G;,v = ^m^" " ^vA^, (4.1) 

for the negatively charged rho meson = pji — ipfi, and the complex conjugate of this equation 
for the positively charged rho meson = + 'P^ ■ Inserting the background gauge field ansatz 

—3 _ 

A2=xieB into the above equation of motion, and Fourier transforming p^ , we obtain the Landau 
energy levels [13] 

(pf T ip2 ) = {eB{2N+l)+pl + mlT 2eB^ (pf ^ ip^ ) (4.2) 

with p3 the momentum of the meson in the X3-direction. The combinations {p^ — /pj") and (pj+ + 
ip^) have spin 5'3 = 1 parallel to the external magnetic field. In the lowest energy state (N = 0, 
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J33 = 0) their effective mass, 

mlff = ml-eB, (4.3) 

can become zero if the magnetic field is strong enough. 

This suggests that the QCD vacuum in a strong magnetic field B = F12 is unstable towards 
condensation of the fields 

p = pf - and p^ = Pj+ + ip:^ 
at a critical value of the magnetic field 

eB, = ml. (4.4) 

Here, in our first approximation, the B-dependence of itself was ignored (as in [2]). AtB^ B^ 
the decay modes of rho mesons into pions, which normally would guarantee a short lifetime of the 
rho mesons, are no longer kinematically allowed because the charged pions' masses increase due 
to the magnetic field (see discussion in [2]). 



mjXGeV-) m^(GeV^) 




eB(GeV^) 



Figure 4: Mass of the rho meson squared as a Figure 5: eBc is determined as the value where 
function of the magnetic field. the effective rho meson mass nt^^j = m^ — eB 

becomes zero, i.e. at the crossing of the 2 curves. 

Next, we can calculate the S-dependence of mp in (4.2) itself, the lowest eigenvalue of the 
equation (2.9), which depends on B both explicitly and implicitly through the changed embedding 
of the probe branes, represented by the value of uq{B). This reflects the effect of chiral magnetic 
catalysis. The result of our numerical study [8] is depicted in Figure 4, showing that mp is an in- 
creasing function of B, as expected. The value of the magnetic field at the onset of the condensation 
becomes shghtly larger: 

eSc = 1.08m^. (4.5) 

We clearly observe that the effect is of the order of a few percent, which is consistent with the 
rough estimate provided in [2]. 

5. Split between deconfinement and chiral symmetry restoration temperature 

In the Sakai-Sugimoto model the chiral symmetry restoration temperature is known to 
increase when a background magnetic field is present [11]. Because of the way the magnetic field 
is introduced in the Sakai-Sugimoto model, the deconfinement temperature Tc is independent^ of 

^Tc is determined solely from the background geometry, which is totally i?-independent in the current approach. 
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B, with Tc 115 MeV for our numerical values of the holographic parameters [8], so there will 
arise a split between the deconfinement transition and chiral symmetry restoration. We numerically 
calculate this split, assuming T^{B = 0) = Tc, and find it to be 1.2% at B = \5m\ k, 0.28 GeV^ and 
3.2% at B = 30m| 0.57 GeV^, see Figure 6. The percent estimation of the magnitude of the 
split is in quantitative agreement with [14, 15], and in qualitative agreement with [10, 16]. Notice 
however that none of these papers work in the chiral limit, and that all of these also witnessed 
an increasing Tc in terms of B. Results in the chiral limit can be extracted from the work [17] 
wherein, interestingly, Tc seems to depend only marginally on B, at least up to eB = 20m^, the 
maximum value considered there. It is perhaps interesting to notice that for ever-growing values 



Figure 6: The chiral symmetry restoration temperature as a function of the magnetic field. The deconfine- 
ment temperature remains at « 0. 1 15 GeV. 

of B, a saturation (horizontal asymptote) is seen in all our quantities {uo{B), mq{B) and T^{B)), in 
agreement with what was found in [1 1]. It is unclear whether this behaviour really corresponds to 
QCD. From QCD-related models hke [10, 14, 16, 18] and lattice simulations like [19], one rather 
expects that e.g. mq{B) (or more precisely, the chiral condensate) and T^{B) will keep growing 
for large B. It looks tempting to associate such a discrepancy to the probe approximation. At large 
values of B, perhaps the quark back reaction might become more prominent, leading to other results 
than found for e.g. niq{B) at large B. Likewise, also the deconfinement temperature might become 
B-dependent beyond the probe approximation, in harmony with works like [10, 14, 16, 18]. It thus 
seems interesting to study the back reaction of the flavour branes on the background geometry, with 
or without magnetic field present^. 

So far, we did not study the actual condensation of the rho mesons, a work which we relegate 
to the future. Let us however already refer to [21] for related work using a different holographic 
setup, and to [22] where the splitting of the 2 transitions was also studied in the Sakai-Sugimoto 
model. However, in contrast to our analysis, no results in physical units were presented, which 
makes explicit comparison with QCD, or QCD models, rather difficult. We did determine all 
"unphysical" parameters of the Sakai-Sugimoto model, like L or Mk, in such a way that known 
QCD physical constants are reproduced as good as possible. It would also be interesting to include 
the pions in a future analysis. Since the rotational symmetry is broken by B = Be^,, one might 
expect a mixing between the longitudinal vector mesons and pions in some manner. It is apparent 
that such happens when taking into account the to the chiral anomaly related Chem-Simons piece 

^In absence of the probe branes, a constant magnetic field B does not change the background geometry, see [20] and 
references therein. 
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of the action [8], but as said before, this effect is suppressed by an additional inverse 't Hooft 
coupling. For sure, getting a clear holographic understanding of the occun^ence of a magnetically 
induced superconducting QCD vacuum at zero temperature, as touched upon in [2, 3], will need 
further efforts. 
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